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Abstract 

We prove several discrete Gagliardo-Nirenberg-Sobolev and Sobolev-Poincare inequali- 
ties for some approximations with arbitrary boundary values on finite volume admissible 
meshes. The keypoint of our approach is to use the continuous embedding of the space 
BV{n) into L'^^'^^-'^\n) for a Lipschitz domain Q C R^, with iV > 2. Finally, we give 
several applications to discrete duality finite volume (DDFV) schemes which are used for 
the approximation of nonlinear and non isotropic elliptic and parabolic problems. 
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1 Introduction 



In this paper, we establish some discrete functional inequalities which are sometimes useful for 
the convergence analysis of finite volume schemes. In the continuous framework, the Gagliardo- 
Nirenberg-Sobolev and Sobolev-Poincare inequalities are fundamental for the analysis of partial 
differential equations. They are a standard tool in existence and regularity theories for solutions. 
The framework is generally used for linear elliptic problems, more precisely it is a classical 
way to prove the coercivity of bilinear forms in i/p , which then allows to apply the Lax-Milgram 
theorem to prove existence of weak solutions. More generally, the framework is crucial for 
the study of nonlinear elliptic or parabolic equations, to obtain some energy estimates which are 
useful to prove existence of weak solutions. Poincare-type inequalities are also one of the step in 
the study of convergence to equilibrium for kinetic equations. 

1.1 Gagliardo-Nirenberg-Sobolev and Sobolev-Poincare inequalities 

In the continuous situation, the Gagliardo-Nirenberg-Sobolev inequality writes as follows. Let 
assume N > 2 and 51 be an open domain of M.^ . Then for 1 < p < and q > I, there exists a 
constant C > such that for all u € W^^PiO.) n L'J{n), 

< ll"lli«(a)' (1) 

where 

P ^11-699 

0<9 < f — < 1 and 



p + q{p—l) m q p N 

We refer to [19l[31] for a proof of this result. We also remind the well-known Sobolev-Poincare 
inequality [TJIS]: 

Wuhnn) < C\\u\\w^^,^nj yu€W''P{n), (2) 

for 

l<q< if 1<P<N, 

N — p 

or 1 < (? < +00 if p > N. 

The mathematical analysis of convergence and error estimates for numerical methods are 
performed using functional analysis tools, such as discrete Sobolev inequalities. Several Poincare- 
Sobolev inequalities have been established for the finite volume schemes as well as for the finite 
element methods. Concerning the finite volume framework, the first estimates were obtained in 
the particular case N = 2, p = q = 2 (which is the standard Poincare inequality) for Dirichlet 
boundary conditions Y. Coudiere and G. Manzini [llj . The idea of the proof in these papers is to 
use some geometrical properties of the mesh. More precisely, given an oriented direction V, any 
cell center of the mesh is connected to an upstream (with respect to V) center of an edge of the 
boundary 9f2 by a straight line of direction T). This connection crosses a certain number of cells 
and their interfaces, and this argument allows to link a norm of the piecewise constant function 
considered with a norm of a discrete version of its gradient. This result was later generalized to 
the case of dimension = 3 by R. Eymard, R. Herbin and Th. Gallouet [16 . Also the same 
method has been applied to get more general Sobolev-Poincare inequalities ^ ioi 1 < p = q <2 
by J. Droniou, Th. Gallouet and R. Herbin 15;, and for p = 2, 1 < g < oo if = 2, 1 < g < 6 
if = 3 by Y. Coudiere Th. Gallouet and R. Herbin [9] , still in the case of Dirichlet boundary 
conditions. Concerning the case of Neumann boundary conditions, a discrete Poincare- Wirtinger 
inequality [p = q = 2) was established in JTHl [20] for A^ = 2 or 3 by using the same method. 



2 



More recently, another idea was used to prove this type of discrete inequahties: the continuous 
embedding of BV{n) into L-^^^^~^^fl) for a Lipschitz domain ft. This argument was first 
exploited in [TB] to prove a discrete Sobolev-Poincare inequality ^ in dimension N = 2 with 
q = 2 and p = 1, in the case of Neumann boundary conditions. Then this method was used in 
[17| to prove general Sobolev-Poincare inequalities ^ in any dimension TV > 1 in the particular 
case of homogeneous Dirichlet boundary conditions. We also mention [S] where the continuous 
embedding of BV{R^) into L^/*^^^^)(R^) is used to establish an improved discrete Gagliardo- 
Nirenberg-Sobolev inequality in the whole space R^, > 1. 

Finally for p — 2, general discrete Sobolev-Poincare inequalities are obtained in [21] for 
Voronoi finite volume approximations in the case of arbitrary boundary conditions by using an 
adaptation of Sobolev's integral representation and the Voronoi property of the mesh. Con- 
cerning the finite element framework, a variant of a Poincare-type inequality {p = q = 2) for 
functions in broken Sobolev spaces was derived in [3] for = 2 and in [71 [35] for = 2, 3. Then 
a generalised result was proposed in [25] , providing bounds on the norms in terms of a broken 
iJi norm (p = 2, 1 < g < oo if iV = 2 and 1 < g < 2N/{N - 2) if iV > 3). The proof is based 
on elliptic regularity results and nonconforming finite element interpolants. Finally, a result in 
non-Hilbertian setting (p ^ 2) was obtained in [13| , taking inspiration from the technique used 
by F. Filbet [TH] and also R. Eymard, R. Herbin and Th. Gallouet fTT , namely the continuous 
embedding of into 



1.2 Aim of the paper and outline 

In this paper our aim is to provide a simple proof to discrete versions of Gagliardo-Nirenberg- 
Sobolev ([!]) and Sobolev-Poincare ^ inequalities for functions coming from finite volume schemes 
with arbitrary boundary values. Several Sobolev-Poincare inequalities are already proved as 
mentioned above but here we propose a unified result. It includes in particular the case of 
mixed boundary conditions. Concerning Gagliardo-Nirenberg-Sobolev inequalities, the result of 
F. Bouchut, R. Eymard and A. Prignet [5] is to our knowledge the only available, and it deals 
with the case of the whole space R^. 

Our starting point to prove these discrete estimates is the continuous embedding of BV{fl) 
into as in [18L 1171 1131 [S] . The main difficulty appears when boundary conditions 

must be taken into account. In the papers mentioned previously [18LI17L[T5] . the boundary condi- 
tions are either homogeneous Dirichlet or Neumann on the whole boundary. In [5 , the problem 
is considered in the whole space R^ . In the case where the function satisfies homogeneous 
Dirichlet boundary conditions only on a part P^ ^ dfl of the boundary, we cannot use the same 
strategy as in [17 , which consists of extending the function considered to R^ by zero. Our idea 
is to thicken the boundary of to take the mixed boundary conditions into account in this case. 

The outline of the paper is as follows. In Section [3 we first define the functional spaces: the 
space of finite volume approximations and the space BV{fl). We will see that BV{il) is a natural 
space to study piecewise constant functions as finite volume approximations. In Section [31 we 
do not take into account any boundary conditions and prove the discrete Gagliardo-Nirenberg- 
Sobolev inequalities (Theorem [3|) and the discrete Sobolev-Poincare inequalities (Theorem \^ 
in this case. These results are the discrete counterpart of ([T]) and ([2]). They may be used for 
instance in the convergence analysis of finite volume schemes in the case with Neumann boundary 
conditions. Then, in Section [H we consider the case where the discrete function is given by a 
finite volume scheme with homogeneous boundary conditions on a part of the boundary. In 
this case, the discrete space (for the finite volume approximations) is unchanged. However, the 
discrete W^'P seminorm will take into account some jumps on the boundary. We prove discrete 
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Gagliardo-Nirenberg-Sobolev inequalities (Theorem[S]) and discrete Sobolev-Poincare inequalities 
(Theorem [5]), similar to ([T]) and ^ but with the V7^'^-seminorm instead of the full M^^'^'-norm. 
Finally, in Section [SJ we show how to extend the results from Sections 12] and S] to finite volume 
approximations coming from discrete dual finite volume (DDFV) schemes. This family of schemes 
is mainly applied to elliptic and parabolic problems. This method can be applied to a wide class 
of 2D meshes (but also 3D |10| ) and inherits the main qualitative properties of the continuous 
problem: monotonicity, coercivity, variational formulation, etc... 



2 Functional spaces 

2.1 The space of finite volume approximations 

We now introduce the discrete settings, including notations and assumptions on the meshes 
and definitions of the discrete norms. Let be an open bounded polyhedral susbset (Lipschitz 
domain) of K^, N > 2, and F its boundary. An admissible mesh of ft is given by a 

family dJl of control volumes, a family £ of relatively open parts of hyperplans in (which 
represent the faces of the control volumes) and a family of points {xK)Kem which satisfy the 
Definition 9.1 in [TB]. It implies that the straight line between two neighboring centers of cells 
{xkjXl) is orthogonal to the face a — K\L. In the set of faces £, we distinguish the interior 
faces a € Smt and the boundary faces a € £ext- For a control volume K g DJl, we denote by £k 
the set of its faces, £int,K the set of its interior faces and £ext,K the set of faces of K included in 
the boundary F. 

In the sequel we denote by d the distance in M^, m the Lebesgue measure in or ^. For 
all <T G £, we define 

d = / '^(^KjXl) for a = K\L e £i„t, 
~ \ d(a:_ftr,(7) for a G £ext,K- 

We assume that the family of meshes considered satisfies the following regularity constraint: 
there exists ^ > such that 

d{xK,cr)>^d„, iorKem, for cr G 5^. (3) 

The size of the mesh is defined by 

h — max (dia.m(K)) . (4) 



In general, finite volume methods lead to the computation of one discrete unknown by con- 
trol volume. The corresponding finite volume approximation is a piecewise constant function. 
Therefore, we define the set X{DJl) of the finite volume approximation: 

X{M) = <u e L^(ri) / 3{uK)Kefm such that = ^ uk^k \ ■ 
[ K&m J 

Let us now define some discrete norms and seminorms on X{dJl). 

Definition 1. Let ft be a bounded polyhedral subset o/R^, DJl an admissible mesh ofil. 
1. For p e [l,-|-oo), the discrete norm is defined by 

Ikllo^p^OT - ( 51 '"Wl^i^in ' yu&x{m). 
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2. In the general case, for p G [l,+oo), the discrete W^'^-seminorm is defined by: 



( 



I U |l,p,95i — 



E 



\a=K\L 



-—J- \UL - UkV 
da 



\ 



,Wu e x{m) 



I 



and the discrete W^'"^ -norm is defined by 

||"||l,p,OT = ||w||o,p,OT + |u|l,p,OT, yueX{M). 



(5) 



3. In the case where homogeneous Dirichlet boundary conditions are underlying (because the 
piecewise constant function comes from a finite volume scheme), we need to take into 
account jumps on the boundary in the discrete W^'^ -seminorm. Let <zT be a part of the 
boundary. In the set of exterior faces £ext, we distinguish f^^^ the set of boundary faces 
included in F". For p G [l,+oo), we define the discrete W^'P -seminorm (which depends on 
TV by 



where 



i.p.ro.OT = 2^ [Dauf , l<p< +00, 
\ae£ J 

\uk-ul\ if cr = K\L e Sint, 
Dau ^ { \uk\ if e £°^t n £k, 

if (7 eSextXS^xt- 



(7) 



We then define the discrete W^'^ norm by 

\\u\\i,p,m = \\u\\o,p,m + |w|i,p,r",OT> Vu e X(97l). (8) 

2.2 The space BV{Q) 

Let us first recall some results concerning functions of bounded variation (we refer to [H for a 
thorough presentation BV{Q)). Let be an open set of and u £ L^{Q). The total variation 
of u in Q, denoted by TVh(u), is defined by 



TVn{u) ^ sup 1^ J u{x)dw{(p{x)) dx, 0eC](f^), < 1, Vx 



(9) 



and the function u G L^{n) belongs to BV{il,) if and only if TVh(u) < +oo. The space BV{n) 
is endowed with the norm 

II u Wsvid) II u Wl^u) + TVn{u). 

The space BV{fl) is a natural space to study finite volume approximations. Indeed, as it is 
proved for instance in [18 , for u E X{dJV), we have 



TVn{u)< J2 

m(cr) \UL 



UK 



|w|l,l,OT < 



<Te£i„t 

(T=K\L 



The discrete space is included in n BV[il). Moreover, ||u||bv(o) < ||m||i,i,ot- 

Our starting point for the discrete functional inequalities is the continuous embedding of 
BV{^) into L^/*^^~^)(ri) for a Lipschitz domain f2, recalled in Theorem[T] 
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Theorem 1. Let be a Lipschitz bounded domain ofM.^, N >2. Then there exists a constant 
c(r2) only depending on f2 such that: 



Af-1 



(10) 



There are also more precise results involving only the seminorm TVn (u) instead of the norm 
II w l|By(n)- Indeed, the seminorm TVh becomes a norm on the space of BV functions vanishing 
on a part of the boundary and also on the space of BV functions with a zero mean value. In 
these cases, the continuous embedding of BV{fl) into L^^^^^^\fl) rewrites as in Theorem [21 

Theorem 2. Let ft be a Lipschitz bounded domain o/R^, N >2. 

1. There exists a constant c(fl) > only depending on such that, for all u G BV{^1), 



n 



\u~u\—dx] < c(f7) TVh(u), (11) 
where u is the mean value of u: 



1 



m{n) 



u{x) dx. 



2. Let C d^l, ^ %. There exists a constant c(r2) > only depending on Q, such that, for 
all u e BV{p?) satisfying u = on T'^ , 



\u\T^dx] <c{n)TVniu). (12) 



Actually, the constant c(ri) involved in Theorems [T] and [2] depends only on 9 and r such that 
the domain has the cone property for these parameters (see [1] Lemma 4-24], [301 Theorem 
8-8]). 

3 Discrete functional inequalities in the general case 

We first consider the general case u G X{^) with the discrete W^'^ norm defined by The 
discrete functional inequalities we will prove may be useful in the convergence analysis of finite 
volume methods for problems with homogeneous Neumann boundary conditions. 

3.1 General discrete Gagliardo-Nirenberg-Sobolev inequality 

We start with the discrete Gagliardo-Nirenberg-Sobolev inequalities which are the discrete coun- 
terpart of ([T]). 

Theorem 3 (General discrete Gagliardo-Nirenberg-Sobolev inequality). Letil be an open bounded 
polyhedral domain ofM.^, N >2. Let 9Jl be an admissible mesh satisfying 

Then for I < p < N and q > 1, there exists a constant C > only depending on p, q, N and Q 
such that 

II u \\o,m.,m < II « \\ip,m II « llor^OT' e ^(^)' (13) 
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where 

o<e< ^ — - < 1 (14) 

- -p + qip-l)- ^ ' 

and 

1 _i-9 6 

m q p N 



Proof. Throughout this proof, C denotes constants which depend only on il, N , p and q. 

As seen in Section [^?^ we have ||'y||5y(Q) < 1,1,50; for aU v G X{9)1). Therefore, applying 
Theorem [U we get 

II V \\o.N/iN-i).m < cm (I V \i.i,m + II v \\o,im) G X{m). (16) 
Let s > 1. For u G X{M), we now define v G X{M) by vk = Iw^l" for all K eM. We note that 



\v\\o,N/{N-i),m= m{K)\uK\"'^' | 

|?^||o,l.OT= Y "^(^)IwkI' = ll"llo..,OT 



IO,s7V/(7V-l),C« 



and 



|w|i,i,OT= X! '7i((t) Im^I" - 



rj=K\L 



But, for all a — K\L, we have 



UK - UL 



Applying a discrete integration by parts and Holder's inequality, we get, for any 1 < p < and 
s > 1: 



a=K\L 



\UK - UL 



Y ^i'^)\uK\'' ^\uk-Ul\ 



K£ma=K\L 



E E ^K--'^ 

Kema=K\L <^ 



E E H<^)da\uK\ " 



(s-l)p 



But, the regularity constraint Q on the mesh ensures that for all K G 9Jl: 



Y < 7 E m(o-)d(a;A',f7) = ^m(A:), 



and then, for any l<p<A'^, s > 1, we get: 



E "^('^^ 



a=K\L 



\uk\ - \Ul\ 



< 



Cs 



(-5-1) 



^{p-l)/p l^ll.P.OT II " ll0,(s-l)p/(p-l),OT 



(17) 
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with 1 < p < and s > 1. Therefore, from we obtain that for u e X{yyi), I < p < N, 
s > 1, and q :— {s — l)p/ (p — 1) > 1: 



u 



However, since 

1 1/s (s-l)/s 

- = — + — < 1, 

s p q 

we obtain by interpolation that 

\\u\\o,sm < Ihllofp.satlhllotg^OT^"' (18) 

which yields 

II 110,^^/(^-1), 9J! 



and finally 

Using the following interpolation inequality 



II II ^ ^ II ll(s-l)/s II l|l/« 

II W||o,sAr/(JV-l),ait S ^(p-i)/(ps) ll"llo,<J,OT ll"lll,p,OT- 



(19) 



lh||o,™,OT < IhllorOiT ll"llJ7an when 1 = + ^ with < a < 1, (20) 

with r = sN/{N - 1), we get 

II u ||o,,„,OT < II w lir,p'OT II u 110^^' VO < a < 1. 

Setting 9 — a/s, with s = (p— l)(7/p+l, we get the expected inequality (jl3p under the conditions 
([Til) and ([13. □ 

3.2 General discrete Sobolev-Poincare inequality 

We now give the discrete counterpart of the Sobolev-Poincare inequalities ([2]). 

Theorem 4 (General discrete Sobolev-Poincare inequality). Let il be an open bounded polyhedral 

domain ofM.^, N >2. Let OK be an admissible mesh satisfying (Q)-([^. 

Then there exists a constant C > only depending on p, q, N and $7 such that: 

pN 

• ifl<p<N,foralll<q<p*:= 



N -p' 



C 

||M||o,q,OT < ^(p„i)/p II "||i,p,aiT, VMeX(OT), 



(21) 



• if P ^ for all 1 < q < +oo 



II u ||o,,,OT < jl^Tjj^ II " lli,P,an, Vu e X(9Jt). (22) 
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Proof. Throughout this proof, C denotes constants which depend only on fi, N, p and q. 
The case p = 1, q ^ p* = N/{N ~ 1) corresponds to estimate Then we obtain (|2ip for 

p — 1 and for all 1 < g < p* by using the fact that (Jl) C L'^(Jl) since fl is bounded. 
Let us first consider the case where 1 < p < N. We get from ^T9\\ that 



\\u\\o,sN/(N-i),m < ^(p_^/(p,) \\u\\o,q,m' \\^\\hp,m (23) 

with 

(s - l)p 

q=- f<l- 

p-1 

sN {s — l)p N - 1 

Then we choose s > 1 such that = , that is s = p, which implies 

A* — 1 p— 1 N — p 

sN pN 



N-1 N-p 

Therefore, we get 

C 

\\u\\o,pN/{N-p)m < ^(p-i)/p \\u\\i,pm yi<p<N, 

and since LP^/(^-p)(f7) C Li{n) for alll < g < pN/{N - p), it yields 

pN 



C 

\\u\\o,qm < ^(p-i)/p \\u\\i,p.m Vge 



1, 



N ~p 



and the proof is complete for 1 < p < N. 

Now for p > N, we begin with the case p — N . Using the Young's inequality in ([T^ . we obtain 



\\u\\n.sN/(N-l),m < C* Q(jv-1)/Ar ll^lll.A'.OT + II llo,(s-l)^V/(N-l),931^ , 



(24) 



with independant of s. Then we proceed by induction on s > — 1 to prove that there exists 
a constant Cg depending on s such that 

C 

II " llo.sAr/(Ar-l),OT < ^(^N-1)/N H " h,N,m- (25) 

For s — N —1, the result is given by ((231) ■ Then let s > such that (pSj) is true for s — 1 . Using 
and we get 



II II ^ C'*(l + Cs_i) 

II " llo,sAr/(JV-l),OT < ^(N-l)/N — ll"lll,W,OT- 

Then ([^5)) is true for all s > iV — 1 , which finally yields 

C 

II M||o,g,OT < ^(Af-i)/Af II "lli,jv,93i V(7e[l,+oo[, (26) 

which is the result for p = N. We emphasize that Cq — > +oo as (7 — > +00. 
Finally for p > TV, we obtain the result using the fact that 

C 

II u \\i,N,m < ^(p_Af)/pAr II " II I'P.g^ yp>N. (27) 
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Gathering ^ and (gT]) we get 

C C 

\\u\\o,q,m < ^(N-i)/N \\^h,N,m < ^^p_i)/p \\u\\i,P,m Vg e [l,+oo[, 

which completes the proof of Theorem 31 

□ 

3.3 Other discrete functional inequalities 

From Theorems [3] and 21 we can deduce a discrete Nash inequahty: 

Corollary 1 (Discrete Nash inequality). Let O be an open bounded polyhedral domain ofM.^. Let 
DJl be an admissible mesh satisfying (3^-^- Then there exists a constant C > only depending 
on ft and N such that 



u 



lo^i ^ ^ II" II 1.2,50^11 ^llo"i,iDP VueX(OT). (28) 



Proof. For N — 2, the result is directly given by the application of Theorem [31 with p = 2, q ^ 1, 
9 = =1/2 and rn = 2. For > 3, let us first apply Holder's inequality: 

II n = Y: ni(i^)|uK|^/(^+^) \u^r/(^^^^ < II u ||;;S+^) II u (29) 
Kem 

Then we apply Theorem H with 1 < p = 2 < and q = p* = 2N/{N - 2): 

C 

II "llo,2W/(W-2),OT < II u||l,2,OT- (30) 

Gathering and dSO]), it yields the result. 

□ 

In the proofs of TheoremjJland TheoremjH we have used the continuous embedding of BV{il) 
into as it is written in Theorem [H But, starting with (|lip instead of (ITUl) leads 

to a discrete Poincare-Wirtinger inequality given in Proposition [T] This result has already been 
proved in [TB]. 

Proposition 1. Let fl be an open bounded polyhedral domain of M.^ . Let 9Jt be an admissible 
mesh satisfying Then there exists a constant C > only depending on fl and N such 

that 

II w||Q,7v/(iv-i),OT < ^^^\u\i,i,m, yuex{m). 
Lt yields estimates only involving the discrete W^'^ -seminorm in some particular cases: 

N 



C 

w-u||o,<j,OT < ;y=|w|l,2,OT e 



N - 1 



\ t \ ^ 

(We recall that u — — -— - / u(x)dx = — tttt > m(K)uK for u G XidJl).) 
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4 Discrete functional inequalities in the case of Dirichlet 
boundary conditions 

In this Section, we consider the case where the finite volume approximation u S is coming 

from a finite volume scheme where homogeneous boundary conditions are prescribed on a part 
of the boundary. This part of the boundary is denoted by r° C F, r° 7^ 0. In this case, the 
natural discrete counterparts of the ly^'^'-seminorm and ly^'P-norm are defined by © and ([H]). 
Moreover, the ly^'P-seminorm becomes a norm on the space of W^'^ functions vanishing on a 
part of the boundary and the Gagliardo-Nirenberg-Sobolev inequalities and the Sobolev-Poincare 
inequalities may be rewritten with the VF-'^'^'-seminorm instead of the ly^'P-norm. Our aim in this 
Section is to prove the discrete counterpart of such inequalities (see Theorem [5] and Theorem [6]). 

As in the general case, the starting point will be the continuous embedding from BV{^) 
into L^/(^~-'^)(n), which rewrites as ([T^ with homogeneous Dirichlet boundary conditions on 
the part of the boundary. However, (fT^ can not be directly applied to u G X(9Jl). Indeed, 
u € X(9Jl) belongs to BV{H) and therefore its trace on the boundary is well defined; but it does 
not necessarily vanish on . Some adaptations must be done in order to apply (jl2p and get its 
discrete counterpart. It will be done in Section l4.1l and yield the discrete functional inequalities 
presented in Section and Section H31 

In this section, we assume the open set is also convex. This will be particularly crucial in 
Lemma [TJ 

4.1 Preliminary Lemma 

We begin with a Lemma which gives the discrete counterpart of (jl2p . This Lemma is crucial to 
prove Theorems [5] and 

Lemma 1. Let Vl be an open convex hounded polyhedral domain ofR^ and ^ ^ be a part of 
the boundary T. Let Wl be an admissible mesh satisfying Then there exists a constant 

c{fl) only depending on such that 

II U llo,JV/(Ar-l),OT < c{n) I U |i,i,ro,OT, Vu e XiM). 



Proof. Let us consider u G since u is piecewise constant, u belongs to BV{fl). Then we 

can define the trace Tu of u by: for almost every x € T, 

1 f 

lim — -— ^ -— / \v — Tv(x)\ dy = 0. 

r^o m(i3(x,r)nf))iB(,,,)na 

Thus in general Tu^-pa ^ and in this framework we cannot take into account the prescribed 
homogeneous Dirichlet boundary conditions u^r = for a E 8^^^ . Therefore the idea is to thicken 
the domain f2 into a larger domain il^ with f2 C fie and define an extension of u to fJg such 
that S BV{U,^) and Tu^ = on a non empty part of the boundary dQ^i which allows to apply 
Theorem [2] to m^. 

Let cr € £ext be a face included in the boundary F. Then cr is a part of an hyperplane % in M^. 
We denote by rio- the unit vector normal to Ti. outward to fi. For every x G M^, there exists a 
unique {y, y') £ H x R such that x — y + y' ricr . For e > 0, we define (see Figure [T]) 

K'^ -.^ {x ^y + y' n„ eR^ : y e cr and <?/'< e} . 
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Since ft is convex, n K^, is empty for all a, a' E Eext with a ^ a' . Now we can define (see 
Figure [1]) 



:= 17 U 



^ U 



The subset Q.^ is polyhedral, then it is a Lipschitz domain. We point out that for all x £ 
A{x, n) < e and then if we consider a face between two new cells and K'^, , we have: 



JV-l 



Then we define 



uY{K%r\Kl,) < Ce 
if 

UK if cr e £ext \ ^exU ^ S Ek-, 



and the function in the following way: 

Ue := ^ ukIk + ^ "i^sl^^S 



We have obviously 

II " llo,A'/(Ar-l),OT < II I1l"/("-i)(0,)- 

Moreover, since the function u^ is piecewise constant and has a finite number of jumps (which 
corresponds to the number of faces a G £ added to the number of faces between the new cells 
K^) we get that belongs to BV{V,s), and 



rV^n.(we) < 5Ini(a)|i^,u|+ ^ Ce 



uk% - Uk-^^ 



Furthermore, since = on a non empty part of the boundary dVl^, we can apply the result 
([T^ of Theorem [21 to Ug. We obtain that there exists a constant c(r2e) such that 



II U ||o,Ar/(Ar-l),OT < c{Vle) 



UK' - UK' 



Now since the open set VL is polyhedral, it is a Lipschitz domain then it satisfies the cone 
condition [1] Definition 4-6] for some cone C and by construction the open set Vl^ also satisfies 
the cone condition with the same cone. Therefore, applying Lemma 4-24 in Jj, the constant 
c(rie) only depends on the dimension of this cone, and not on e > 0. Then passing to the limit 
e — )■ we finally get that 

\\u\\o^N/(N-l)m < c(ri) I U |i,i^OT. 

□ 

Now using this Lemma we can prove the discrete Gagliardo-Nirenberg-Sobolev and Sobolev- 
Poincare inequalities in the case with some homogeneous Dirichlct boundary conditions. 
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Figure 1: Construction of the cell (left) and of the domain £7^ (right). 



4.2 Discrete Gagliardo-Nirenberg-Sobolev inequalities 

Theorem 5 (Discrete Gagliardo-Nirenberg-Sobolev inequalities). Let U, be an open convex 
hounded polyhedral domain of and ^ % he a part of the houndary. Let dJl he an ad- 
missihle mesh satisfying Then for any I < p < N and q > 1, there exists a constant 

C > only depending on p, q, N and such that 

C 

II u ||o,™,OT < ^(p-l)0/p I w li,p,ro,OT II w \\l~q,m^ e X{mi), (31) 
where 9 and m satisfy (|14p and (jl5p . 

Proof. The proof is similar to the proof of Theorem[31 Let 1 < p < N and s > 1. For u e X(9JT), 
we apply Lemma [T] to v € X(9Jt) defined by vk = u|f for all K g 971. It yields 

ll"llasAf/(Ar-l),OT - ^(p_i)/p l'"kp,r",OT II llo,(s-l)p/(p-l),93t 

with 1 <p < N and s > 1. 

Choosing q = (s — l)p/ (p — 1) (s = 1 + g — - > 1), we obtain 

P 



C 

U\\o^sN/{N-l),<m S ^(p_i)/p^ l"li;p,ro.OTll"lld,9,OT 



< 77;;^l"lliT0,..zll"llL:^^^ (32) 



Then, using an interpolation inequality as in the proof of Theorem [3J we get 

II II ^ II 11'^ II iii^'^ ^ ^ I 1^/^ II \\^~{^/^) 

||U||0,m,OT S II " llo,sAr/(Ar-l),OT II " llo,q,OT - ^a(p-l)/ps l'"ll,p,ro,OT ll"llo,g,OT ' 

Taking B — aj s concludes the proof. 



□ 
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4.3 Discrete Sobolev-Poincare and Nash inequalities 

In the case with some homogeneous Dirichlet boundary conditions, the discrete Sobolev-Poincare 
inequahties rewrite as follows. 

Theorem 6 (Discrete Sobolev-Poincare inequality). LetQ be an open convex bounded polyhedral 
domain ofM.^. Let 971 be an admissible mesh satisfying (0^-^. Let C P, with P*^ ^ 0. 
Then there exists a constant C > which only depends on p, q, N and such that: 

pN 

• ifl<p<N,foralll<q<p*- 



N ~p 
C 



c 

u|lo,g,OT < 77— Tw- |w|i,p,ro,OT yueX{m), (33) 



j/p ^ N, for all 1 < q < +oo 



c 

\\u\\o,q,m < c(p-i)/p l^li.P,r»,9Ji yueX{m). (34) 



Proof. The proof is similar to the proof of Theorem HI starting from (|32l) instead of . □ 

Now using Theorems [5] and [51 we easily get a discrete version of Nash inequality: 

Corollary 2 (Discrete Nash inequality). Let be an open convex bounded polyhedral domain of 
M^. Let DJl be an admissible mesh satisfying Let C P, with P'^ 7^ 0. 

Then there exists a constant C > only depending on 17 and N such that 

II " lloji ^ I " II " llo*i,an e (35) 



5 Application to finite volume approximations coming from 
DDFV schemes 

The discrete duality finite volume methods have been developed for ten years for the approxi- 
mation of anisotropic elliptic problems on almost general meshes in 2D and 3D. They are based 
on some discrete operators (divergence and gradient), satisfying a discrete Green formula (the 
"discrete duality"). The DDFV approximations were first proposed for the discretization of 
anisotropic and/or nonlinear diffusion problems on rather general meshes. We refer to the pio- 
neer work of F. Hermeline [221 [23 IMl 123 IH] who proposed a new approach dealing with primal 
and dual meshes and Y. Coudiere, J. -P. Vila anf Ph. Villedieu [11 who proposed a method of 
reconstruction for the discrete gradients. Next, K. Domelevo and P. Omnes [14], S. Delcourte, 
K. Domelevo and P. Omnes [T^ presented the discrete duality finite volume approach (DDFV) 
for the Laplace operator. Then, B. Andreianov, F. Boyer and F. Hubert |2j gave a general back- 
ground of DDFV methods for anisotropic and nonlinear elliptic problems. Most of these works 
treat 2D linear anisotropic, heterogeneous diffusion problems, while the case of discontinuous 
diffusion operators have been treated later by F. Boyer and F. Hubert in [6]. F. Hermeline 
[23 [2S] treats the analogous 3D problems, S. Krell [2Z] treats the Stokes problem in 2D and in 
3D whereas Y. Coudiere and G. Manzini [TIP treat linear elliptic convection-diffusion equations. 
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The construction of DDFV schemes needs the definition of three meshes: a primal mesh, a 
dual mesh and a diamond mesh. Then, the approximate solutions are defined both on the pri- 
mal and the dual meshes, while the approximate gradients are defined on the diamond meshes. 
Therefore, we need to adapt the definition of the spaces of approximate solutions and the defi- 
nition of the discrete norms. It will be done in Section 15.11 Then, we will be able to establish 
some discrete Gagliardo-Nirenberg-Sobolev and Sobolev-Poincare inequalities, in the general case 
(Section 15. 2p as in the case with Dirichlet boundary conditions (Section 15. 3p . 



5.1 Meshes and functional spaces 

Meshes. Let fl be an open bounded polygonal domain of K^. The mesh construction starts 
with the partition of fl with disjoint open polygonal control volumes. This partition, denoted by 
971, is called the interior primal mesh. We then denote by ddJl the set of boundary edges, which 
are considered as degenerate control volumes. Then, the primal mesh is defined by 071 = 97tU997l. 
To each primal cell K. S 97t, we associate a point x/c G IC, called the center of the primal cell. 
Notice that for a degenerate control volume /C, the point x/c is necessarily the midpoint of IC. 
This family of centers is denoted by A" = {xac, /C G 971} . 




Interior and exterior 
primal meshes 



Interior and exterior 
dual meshes 





V J \ ■ /! 


'/ ' 





Diamond mesh 



Figure 2: Presentation of the meshes 

Let X* denote the set of the vertices of the primal control volumes in 97t. Distinguishing 
the interior vertices from the vertices lying on the boundary, we split X* into X* = X*^^ U X*^^. 
To any point x/c* € '^hiti associate the polygon /C* obtained by joining the centers of the 
primal cells whose X)c* is a vertex. The set of such polygons defines the interior dual mesh 
denoted by 971*. To any point x^c G X*^^, we then associate the polygon /C*, whose vertices are 
{xk.*}^{xk. G X /x]c* e e 97t}. It defines the boundary dual mesh 9971* and the dual mesh 
is defined by M* = 971* U 9971*. 

In the sequel, we will assume that each primal cell /C G 971 is star-shaped with respect to xx: 
and each dual cell /C* G 971* is star-shaped with respect to Xjc* . 

For all neighboring primal cells K, and iZ, we assume that dJCndC is a segment, corresponding 
to an edge of the mesh 971, denoted by cr = JC\C. Let £ be the set of such edges. We similarly 
define the edges £* of the dual mesh 971*: a* = IC*\C*. For each couple (ct, cr*) e £ x £* such 
that a = )C\C = ixic*,xc*) and cr* = IC*\C* ~ (x/cxc), we define the quadrilateral diamond cell 
Vcr^a" whose diagonals are a and a* . If ct G £ n dH., we note that the diamond degenerates into a 
triangle. The set of the diamond cells defines a partition of fl, which is called the diamond mesh 
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and is denoted by T>. Let us note that 33 can be splitted into S) = 'Dint U T)ext where 'Dint is the 
set of interior (non degenerate) diamond cells and Dext is the set of degenerate diamond cells. 

Finally, the DDFV mesh is made of the triple T — (9Jl, dJl*,D). See Figure[2]for an example 
of DDFV mesh. 



\ 




> \. 



/ 

■ Vertices of the primal mesh \ 

* 

o Centers of the primal mesh \^ 

— a = IC\C, edge of the primal mesh \ 

* 

— (7* = IC*\C* , edge of the dual mesh \ 
Diamond 2?„ „■* 



Xc 



oxc :: 



dK',i 



lex 



Figure 3: Definition of the diamonds "D^jy* 

Let us now introduce some notations associated to the mesh T. For each primal cell or 
dual cell V in or 9Jt*, we define my, the measure of the ceU V, £v^ the set of edges of V, 
Dv = {2?a,o- e S), cr e dy, the diameter of V. For a diamond I?ct,(t*, whose vertices are 
{XK ), we define : m^ and mo-* the lengths of the primal edge a and the dual edge 

a*, Taj), the measure of V, dv its diameter and the angle between {xic,xc) and (a^K* , a;£-» ). 
As shown on Figure [31 we will also use two direct basis (tk;*^£», iIo-k) and (no-'c* j ''"tc,^), where 
iIctk; is the unit normal to cr, outward /C, rio-'K* is the unit normal to a*, outward /C*, Tk.*,ii* 
is the unit tangent vector to u, oriented from /C* to £*, Tk.,c is the unit tangent vector to u*, 
oriented from K, to C. For a boundary edge a = [xfc , xc] G we define djc*^c the length of 
the segment [xic',xc] and dc*,c the length of the segment [xc'^xc]- 

In all the sequel, we will assume that the diamonds cannot be flat. It means : 



3ar €]0, — ] such that | sin(ax>)| > sin(ar) V2? G S). 
As for all 1)0,0* G D, we have 2mu = mcrmcr* sin(ax)), the hypothesis p6p implies 



m„m^. < 



2mx) 



sin(Q;7-) 

We also assume some regularity of the mesh, as in ^ , which implies 



(36) 



(37) 



3C > 0, nvm^. < ^ V/C e mt, 



c 



E 



m^m„. < 



(38) 
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Definition of the approximate solution. A discrete duality finite volume scheme leads 
to the computation of discrete unknowns on the primal and the dual meshes : iu)c)ic£M ^^'^ 
{''J-K')ic-'£'m'' From these discrete unknowns, we can reconstruct two different approximate solu- 
tions : 

But, in order to use simultaneously the discrete unknowns computed defined on the primal and 
the dual meshes, we prefer to define the approximate solution as 



Therefore, the space of approximate solutions Z{T) is defined by: 



such that u 



2 

For a given function u G Z{T), we define the discrete norm by 




Discrete gradient. A key point in the construction of the DDFV schemes is the definition 
of the discrete operators (divergence and gradient). We just focus here on the definition of the 
discrete gradient, which will be useful for the definition of the discrete W^^'^'-seminorms. 

Let u S Z{T). The discrete gradient of u, V^u is defined as a piecewise constant function on 
each diamond cell : 

•Pes 



V w = r I n^rK H no 



where, for I? g £>, 

sin(ai,) 

This discrete gradient has been introduced in 11 . It verifies: 

V u-Tk',i:' = and V u-Tk,c. = ■ 

mo- mo- 
Using this discrete gradient, we may now define the discrete VK^'^-seminorm and norm of a given 
function u £ ZiJ'): 



i/p 

l"li:P,r = ( 

||u||i,p,r = ||u||o.p,r + |u|i,p,r 
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5.2 Discrete functional inequalities in the general case 



Our aim is now to extend the results of Section [3] to the case of finite volume approximations 
coming from some DDFV schemes: u £ Z{T)- We will use that such functions are defined as 
u — -^(uot + Uppp) with Wot G X(9Jl) and G X(9Jt*). Nevertheless, we may take care because 
the primal and the dual meshes St and 931* does not satisfy the admissibility condition required 
in Theorems |3] and H) 

Theorem 7 (General discrete Gagliardo-Nirenberg-Sobolev inequality in the DDFV framework). 
Let Q, he an open hounded polyhedral domain of M.^ . Let T = (9Jl, 9Jt*, S) he a DDFV mesh 
satisfying ([55)1 and (155)) . 

Then for \ < p <2 and g > 1, there exists a constant C > only depending on p, q and such 
that 

""""'-^ - (sin(a^))^/^C^(^-^)/(^) ll"ll?,,,rll"lli:Jr. yuezin (39) 

where 

< 6» < ^ < 1 

p + q[p- 1) 

and 

11-0 6 
m q p 2 



Proof. We start as in the proof of Theorem[21 Let s > 1 and set q — {s — l)p/{p — 1) > 1. For 
u e Z{T), as u<xn € X(93t) and G X{DJl*), we may write: 



0,2s, an 



< c( 



um\ 



li,i,an 



U: 



But, following the same computations as in the proof of Theorem [31 we get 



1 1,1, OT 



J2 



\uk\ ~ \uc\ 



< s ^ mcrmcr* 



UK - Uc 



m„ 



{\uKr' + \ucr'] 



(40) 
(41) 



< s 



E 



m„ 



Using the regularity hypotheses on the mesh, we get 



usm\ 



C 



UK ~ Uc 



E 



E E m^m^HwK;! "-^"^ 



m-p 





• 




nicr* 





But, by definition, — — — V^u • Tk.,c and therefore 



m^* 



UK - Uc 



m„ 



< |V^u|. It yields : 



C 



"^'"i^i'^- (sin(ar))i/PC(^'-i)/^''"'''^''^""'""°'''™' 
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Injecting this last inequality in (j40p . and using the interpolation inequality (jl8p . we get: 

C 

\M\l2s,m < (sin(a^))i/pC(p-i)/p ll^^^llorq!^ (I^kp.r + II"otIIo,p,ot) 

C 1 

because ||wotI1o,p,ot - ^ II 

^lli p T definition. Then, using the interpolation inequality (j20p with 
r = 2s and — = ^—^ + — with < a < 1, as in the proof of Theorem [31 we obtain: 



C 



r ii"llo,g,r 



With similar computations on the dual mesh, from (PT|) . we get 

ll"SFllo.m,SF - (sin(ar))"/("'')C"(P"^^/("P) H^H lir ll"llo,g,r • 

Finally, setting 9 = a/s with s — {p — l)q/p + 1, it yields the expected inequality (15^1) . □ 

As in the classical finite volume framework, we can now prove discrete Sobolev-Poincare 
inequalities. The proof is similar to the proof of Theorem |4l it will not be detailed here. 

Theorem 8 (General discrete Sobolev-Poincare inequality in the DDFV framework). Let il be 
an open bounded polyhedral domain o/K^. Let T = (9Jl, D) be a DDFV mesh satisfying pB]) 
and ([55)1 . 

Then there exists a constant C > only depending on p, q and il such that: 

2p 



• if'^<P<2,foralll<q< 



2-p' 



^ 2, for all I < q < +co, 



C 

(sin(a^))VPC(P-^)/W V^eZ(r). (43) 



Let us now focus on the Poincare-Wirtinger inequality in the DDFV case. This result has 
been proved recently in ^29,. We will give here a proof using the embedding of BV{fl) into L'^{il) 
dTTl) recalled in Theorem!^ 

Theorem 9 (Discrete Poincare-Wirtinger inequality in the DDFV framework). Let ft be an 
open bounded polyhedral domain ofM.'^. Let T — 971*, S) be a DDFV mesh satisfying (|36p. 
There exists a constant C > depending only on fl, such that for all u G Z{T) satisfying 

m^UK; = ^ ni^.u^- =0, (44) 

we have 

ll-llo,..r<-^Hi,2,r- (45) 
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Proof. Let u e Z{T). Applying ((TTl) to e X{M) and G we get, under the 

hypothesis 

ll"OTllo,2,Et)i ^ c(f7)rVs2(uOT) < c{n)\u<M\iirgi 



But, 



II ^ \UK-Uc\ 



sin(Q;7-) 



1,2,T ' 



thanks to Cauchy-Schwarz incquaUty. By the same way, we get the same bound for jurg^lj^ ^ 

2 

and it finally yields ||u|jQ 2 7- — ym(ri)^/^c(f7) 2 r- '-' 

5.3 Discrete functional inequalities in the case with Dirichlet boundary 
conditions 

In this Section, we want to extend the discrete Gagliardo-Nirenberg-Sobolev inequalities of Sec- 
tion 14.21 to finite volume approximations obtained from a DDFV scheme. We first recall how 
Dirichlet boundary conditions are taken into account in DDFV methods. Let Fq be a non empty 
part of the boundary. At the discrete level, homogeneous Dirichlet boundary conditions on Fq 
will be written: 

uic = 0, V/C € dm, JCCT° and ujc- = 0, V/C* G dm*, WdT" ^(d. (46) 

Therefore, we consider the corresponding set of finite volume approximations, Z^o^T) defined by: 

^ro(T) ={ueZ{T) satisfying gSI} ■ 

Let us note that the definition of the discrete W^'P seminorm is the same on Zro(T) as on Z{T)- 
Indeed, the fact that the approximate solution vanishes at the boundary is taken into account 
in the definition of the discrete gradient V^w for V g 2)ext, and therefore in \u\-^ ^ ^. 

Finally, combining the techniques of proof of Theorem [5] (using Lemma [Ij and Theorem [7j 
we establish the following Theorem. 

Theorem 10 (Discrete Gagliardo-Nirenberg-Sobolev inequality in the DDFV framework). Let 

57 be an open convex bounded polyhedral domain of and F" be a part of the boundary. Let 
T = (M,M*,£>) be a DDFV mesh satisfying ^ and 

Then for 1 < p < 2 and q> 1, there exists a constant C > only depending on p, q and fl such 
that 

(J 

ll^llo.™^r < (sin(a^))g/PC^(p-i)/(p) l^lLrll^lltg'r' VueZro(r), (47) 



where 

(X < 

p + q(p-l 



O<0< ^ < 1 
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